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Abstract 

A solution with the pole conhguration in six dimensions is analyzed both 
analytically and numerically. It is a dimensional reduction model of Radall- 
Sundrum type. The soliton configuration is induced by the bulk Higgs mech- 
anism. The boundary condition is systematically solved up to the 6th order. 
The Riemann curvature is finite everywhere. An exact solution for the no 
potential case is also presented. 
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1 Introduction 



Seeking a realistic model of our world, from the higher-dimensional stand- 
point, the brane- world approach is now vigorously envestigated in various 
ways. It started with the wall configuration in 5 dimensions(D) space- 
time |]l|, |[. Basically the brane- world physics makes use of the soliton config- 
uration. We can explore new physics, including the beyond-standard model, 
by making use of the distinguished properties of soliton, such as the local- 
ized zero modes, stable configuration, singularity-free behavior due to the 
extendedness, etc.. From the reason of improving some difficulties in 5D 
model[^, ^, ^, or of getting more general aspects of the brane- world physics, 
further-higher-dimensional models ( string-like object in 6D, monopole con- 
figuration in 7D, • • •) are also envestigated]^, |^, |8|, ^, |10|. In this paper, we 
examine a 6D model. 

In the original model, 3-brane(s) is introduced, by hand, as a 5-function 
distribution in the extra space. We do not take such approach in order to 
seek a non-singular solution. We induce the same configuration as a soliton 
(kink) solution. This situation is very similar to the relation between the 
Dirac and 'tHooft-Polyakov monopoles. The latter is the soliton (in the 
gauge+Higgs system) interpretation of the former system where the gauge 
potential has 5-function distribution in the plane perpendicular to the Dirac 
string. We seek a soliton solution in the 6D gravity + Higgs system. The 



main purpose is to establish the pole solution |11|. The most crucial point is 



to confirm the convergence of some infinite series appearing in the solution. 



It guarantees the boundary condition. In [l^] they are solved up to the 2nd 
order approximation. Here we show them up to the 6th order. The result 
strongly shows that the series converge and the boundary conditions are 
sufficiently satisfied. 

2 The Dimensional Reduction Model in Six Di- 
mensions 



In 1 11 1 a solution with the pole configuration in six dimensions is presented. 



The model is the 6D gravity coupled with the Higgs fields. 
S[Gab,'^] = J d'^XV^i-^M^R - G^^dA<^*dB'^ - V{^*, $)) 



y($*,cI>) = ^(|cl>|2_V)2 + A , (1) 
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where (X^) = {x^, p, ^p), ^ = 0, 1, 2, 3. x^'s are regarded as our world coor- 
dinates, whereas {p, if) the extra ones( < p < oo, < (p < 2tt). Gab is the 
6D metric field, <I>* and $ are the complex scalars (Higgs fields). M{> 0) is 
the 6D Planck mass. A(> 0),vo{> 0) and A in the potential V are vacuum 
parameters. Expecting Randall-Sundrum type dimensional reduction, we 
assume the following form as the line element. 

ds^ = e-^'^'^P'^rif^^dx^dx" + V + p2e-2-(p)^^2 ^ 

< p < oo , < 99 < 27r , (2) 

where r/^jy = diag(— 1, 1, 1, 1). This is a natural 6D minimal-extension of the 
original 5D model by Randall-Sundrum Q. Two "warp" factors, e"^"^*-''^ 
and e~'^'^^p\ appear. For the fixed p case {dp = 0), the metric is the Weyl 
transformation of the product-space, the 4D Minkowski x the circle S^. 
The coordinate p can be regarded as the scaling parameter. We require 
the periodicity in ip for ^{p,ip}) and assume the form: $ = P(/?)e*™'^ (m = 
0, ±1, ±2, • • •). The 6D classical field equations of (|l|) are obtained as 

3^' _ ^ + a'u;' + Lo" - (uj'f + 2-= 2M~^P'^ , 
P P 

-16cj'^ 4— - Aa'J + 4a" = 2M~^V , 
P 

ia" -Wa'^ = M-\P'^ -m^^P^ + V) . (3) 

p^ 

We take a unit M = 1 in the following for simplicity. ^ 



3 An Exact Solution for the No Potential Case 

In order to see the structure of the solution, it is useful to consider a simple 
case, that is, the no potential case: A = A = {V{^*,^) = 0). Even in 
this case, the scalars are still coupled with gravity. In this section we stress 
those properties which are shared with the general case of Sec. 4. We consider 
the m = (no flux) solution. The following exact solution of the rational 
function type can be obtained. 

// N 1 ^ + 5 ^, ,,2 2(2A-5) 



^ This is allowed by the following invariance of (0) or , due to the simple dimensional 
counting: $ ^ M^$, (P ^ M^P) , vo ^ mKq, A M'^X ,A M'^A , 
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where A{> 5/2) and B are integration constants. j4 is a dimensionless 
constant, whereas B is a dimensional one ( the dimension of length ). The 
lower bound for A comes from the positivity of P'^. This is a 2-parameters 
family solution (for a' ,uj' , and P') at this stage. Some special cases are i) A = 
5/2, P' = 0; ii) B = 0, -\a' = J = {4/A)xp-\P'^ = {2{2A-5)/A^)x p-^; 
iii) A = i, J = p^^\ iv) A ^ oo, ct' = 0, = 0,P' = 0. The integration 
constants are fixed by a required boundary condition. 

Due to the condition A > 5/2, we do not find the solution which satisfies 
the RS-type configuration: a' const. (> 0) as p — > oo. However, if we 
restrict the region of p as \B\ ^ Ap, then we have a' ~ —1/B, to' ~ 
(l/p) - {A - A)/B, P'^ ~ 2(2^ - 5)/B^. This looks a RS-type metric 
although we cannot take p ^ oo. 

The final solution for a, uj and P are given as 

1 \Ap + B\ \p\.\Ap + B\-^ 

^(■p) = -T 1^ — f; — ' ^(^) = 1" ■ 



where C(> 0),D{> 0),E{> 0) are another integration constants with the 
dimension of length. They come from the symmetry of the constant transla- 
tion {a — > a" + const.,a; — > a; + const., P — > P + const.) in eq.(^) with V = 0. 
This final form (^) is a 5-parameters family of solutions. 
The line element is given by 

^ + + l^i^f-'"' . V . 

2 4 6 2M-4) 

For the case of ^ ^ oo the space-time becomes locally 6D Minkowski with 
the deficit angle 6 = 27r(l — A) = — oo. For the restricted region \B\ S> Ap, 
the above metric reduces to 

ds^ = e s P I ^ I rt^^dx^'dx" + + | ^ | • i^^^^a _ 

As a form, this is a RS-type metric although we cannot keep it in p — > cxd 
region. 

The 6D Riemann scalar curvature is obtained as 
^--{Ap + Bf-^ ' 
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which is negative semi-definite. The absolute value goes to as p — > oo. For 
A = 5/2 (the case i) above), R = is valid everywhere. The metric in this 



case 



2 ' 

r|5„, 0,^4/5 ( >|6/5 

^,2 ^ I \2P+^ I ^^^dxV + dp' + l I • E'd^' . (9) 

For oo > A > 5/2, we can see three cases: a) i? > 0, There is no curvature 
singularity; h) B = 0, The curvature is singular at p = 0; c) S < 0, The 
curvature singular at p = —B/A. As for the horizon, from the condition 
< ;! < |, there is no horizon in the region < p < oo. Hence, in b) and 
c), the singularities are naked ones. 

The structure of the solution (|6l) is similar to the Kasner solution [1121 in 



1+3 D. Ref.[13| explains the cosmic string solutions in the Abelian Higgs 
model in 1+3 D. The two branches explained there, that is, the cosmic string 
branch and the Melvin branch look to have correspondence with the A ^ oo 
case and A = 5/2 case respectively. 



4 Pole Solution in Six Dimensions 



Now we consider for the case of the spontaneous breakdown (Higgs) poten- 
tial: A > 0,^0 > 0, and A is general(at this stage). The boundary condition 
for P should be taken as 

p^oo , P{p)^+vo , (10) 

where we take the plus sign, using the freedom of $ ^ — ^> in (^, for the 
asymptotic value of P. We require (based on the analogy to 5D Randall- 
Sundrum model [I], ^, H) that a' and go to constants as p — > oo. Then 
we can deduce, using (y), that 

m = (no flux) , 
p ^ oo , a' ^ a , Lo' ^ a , a = +^ ^^M^^ . (11) 

This result says A should be taken negative: A < 0. We choose the plus 
sign as the asymptotic value so that the Weyl scaling factors in (^ work 
as shrinking the proper distance as p increases. This choice guarantees the 



5 



stableness of the system. As for the behavior near p = (ultra-violet region), 
we can take, in a consistent way with (|3|), as 

/9 ^ +0 , a' ^ sp" , J ^ wp^ , P ^ xp" 
a = 5 = l, c = , s/0,u;7^0,x/0 , 

s = ^^{x^-vo'? + \ , (12) 

with the additional possible = with c > 1. 

We examine the following form of solutions 1 11]. 

oo 

^'(P) = « E 7^^{t-nh(M}^"+^ , 

oo T 

-'(^) = «E(|rTI)!^*""'^^^^)^'"^' ' 

oo 

^(z') = ^0 E TSjita^MMP" , (13) 

n=0 ^ ' 

where c's, d's and e's are coefficient-constants to be determined, a' and u)' 
are composed of odd powers of tanh(A;/9), whereas P is of even powers. They 
come from the behavior at the ultra-violet region ([l^). A; is a new scale 
parameter which shows the thickness of the pole. In the following, we take 
k = 1 for simplicity. ^ 

The key equation for obtaining the solution of the form ([l^) is the fol- 
lowing expansion formula about 1/ p. ( The factor 1/ p appears in @ ) . 



= S2n , 

x=0 



1 2 °° S2n d"^^ ( X 

~p " (2^^*''''^ ' d^ vi^i?f 

1 1 12 

< p < oo , So = - , S2 = -- , S4 = — - , ••• , (14) 

2 3 5 

All coefficients are finite as shown above. We can take the limit p ^ oo 
above, and see the infinite series converges and gives : J2'?^=o (^)T ~ 0. We 
expect the coefficient-series appearing in (13) have the similar behaviors. 
It is a key, in the present treatment of the infinite series, that we take the 
expansion using powers of tanh p not those of e"^''. Q 

^ This is because eq.(^ is invariant under the change: p kp, X — > (l/fc^)A, A 
(l/fc^)A, vo^vo. 

* We can not reexpress the RHS of the first equation of ( |l4| ) as 1/p — X/^o (2n)! 
with finite coefficients. The same notice is said about (p^. 
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We can show the above form of solutions indeed satisfy (^) if the coef- 
ficients satisfy some recursion relation 1 11]. The first few orders are given 
as 



C3 



ci = i^{lXvoHl-eoY + M 
di = -|ci , 

Co : free parameter, 

3 X^VQ^ „ 2/1 ^ 2\2 



32^^60^(1 - + ci(2 + 5aci) 
ds = -|ci(l + baci) , 
62 = -|At;o2eo(l -eo^) • 



(15) 



We note one free parameter cq appears [| and all coefficients are expressed 
by 4 parameters {X,vo, A,eo). The parameters, however, have 3 constraints. 



n=0 



C2ra+1 
(2n + 1) 



J^2n+1 



E e2n 
(2n)! 



o(2n + l)! 



(16) 



n=0 



which comes from the boundary conditions ( |To| ) and ([Tl|). Therefore the 
present solution is a l(=4-3) parameter family of solutions. 



5 Evaluation of Coefficients and Analytic Result 

We present here the results of concrete values of c's, d's and e's for an input 
value eo(= P{0)/vo) = —0.8. We solve constraints (|l6|) by taking the first 7 
terms (up to n=6th order). The most crucial point of the present solutions is 
to confirm the convergence of the infinite series ^^"^=0 (2n+i)! ' Sn=o (tn+iy. ' 
and J2'^=o ^^Jyi' ■^hich guarantees the present boundary condition. The 6th- 
order approximation calculation gives the vacuum parameters as 

vo = 0.9625 , A = -2.725 , A = 18.375 . (17) 

f\ f\ The obtained values of the coefficients are shown in Fig. 1 . In the figure 

^ This freedom corresponds to, when A = 0, the constant translation symmetry (P 
P + const.) in eq.(|^ . 

® These values should not be confused with numerical results like Sec. 6. They should, 
in principle, be definitely determined just like the energy eigenvalues of the hydrogen atom 
are fixed by the boundary condition. 

^ The constraint criterion is (1 - ^Lo + (1 " ELo + (1 " 
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Fig.l The values of C2„+i/(2n + l)!(large blob),(i2n+i/(2n + l)!(small blob) 
,e2n/(2?^)! (small circle). The large circles show (0.8)'^. ( n = 0, 1, • • • , 6 .) 



we also plot the data from the geometrical series: 1/(1 — x) = 1 + x + ■ ■ ■ 
at X = 0.8. Comparing them, we can recognize the convergence of the 
coefficient-series (up to this approximation order). Note that all three series 
are 'oscillating', which is advantageous for their convergence. Using these 
results, the analytical results of P{p),a'{p) and io'{p) ( (|l3|))are shown in 
Fig.2. We can see comparative behaviors of two 'warp' factors {a',uj') in the 
ultra violet region, p < 2. In particular we notice a very 'delicate' dip in 
the behavior of the second 'warp' factor u' sX kp ^ 1. The (6D) Riemann 
scalar curvature is also shown in Fig. 3. It shows the curvature is positive 
both inside and outside of the pole, whereas negative (small absolute value) 
or zero between them. It is non-singular everywhere. There is no horizon in 
< p < oo. The deficit angle at p = c» is 27r. 



6 Numerical Results 

As the coupled differential equations for P{p), cr'{p) and uj'{p), the equations 
(^) with m = 0, have the standard form of the numerical analysis, that is, 
Runge-Kutta method. We can numerically solve (^ without any ansatz 
about the form of the solution. In this approach, the choice of the initial 
values with high precision is required. In the present case we cannot take 
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Fig. 2 The analytic result of u;'/a(bold line), cj'/a(dashed line) and 
P/t^o (normal line). The graphs are depicted by using (ll) with the 6-th 
order approximation. The horizontal axis is kp. We take k = 1. 
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k rho 



Fig. 4 The numerical results for u;'/a(top), cr' /a (middle) and P/?;o(down). 
They are obtained by Runge-Kutta method. One step value along fcp-axis 
is 0.025. About 65 points are plotted for each line in the figure. The initial 
point is p = 0.1. The horizontal axis is kp. We take k = 1. 



p = as the initial point because the appearance of the factor 1/p in (^). 
(Note this does not say the solution is singular at p = 0. We seek a non- 
singular solution.) We take /O = 0.1 as the initial point of p. As for the initial 
values we borrow the results from Sec. 5 : P(O.l) = —0.79386, (t'(O.I) = 
-0.10488, uj'{Q.l) = 0.070046. The result by the Runge-Kutta calculation 
is shown in Fig. 4. It shows the analytic solution of Sec. 5 is well reproduced 
in the ultraviolet region ( small-/? region ) but not for the infrared region. 
We should notice that the the consistent region (with the analytic result) , 
p < 1.3, becomes definitely larger than that of the 2nd order approximation 
[11 1 , p < 0.9. The numerical output data stop at p ~ 1.7 with producing 
imaginary values. This occurs because keeping the positivity, P'{p)^ > 0, 
becomes so stringent in the infrared region. The quantity becomes so small 
in that region and vanishes at /O = oo. We understand that further higher 
precision is required for the initial values in order to extend the consistent 
region furthermore. 



7 Conclusion 

We have presented the both analytical and numerical solution in the 6D 
reduction model. The approximation oder is n = 6. Compared with the 
previous result of n = 2 order one, the results improve definitely. Much 
evidence about the convergence of the coefficients series is obtained. An 
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exact solution, for the no potential case, is also obtained. 

We notice rather stable behavior of the following values under the in- 
crease of the approximation order : the vacuum parameters (fo,A,A); the 
coefficients (c2n+i, d2n+i, e2n\'n = 0, 1, • • • , 6); the form of a' , uj' and P. This 
strongly indicates the result of Fig. 2 is quite near the exact result. 
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